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A new cosmological variable is introduced which characterizes the degree of departure from Ein- 
stein's General Relativity (GR) with a cosmological constant. The new parameter, zu, is the cos- 
mological analog of 7, the parametrized post-Newtonian variable which measures the amount of 
spacetime curvature per unit mass. In the cosmological context, w measures the difference between 
the Newtonian and longitudinal potentials in response to the same matter sources, as occurs in 
certain scalar-tensor theories of gravity. Equivalently, zu measures the scalar shear fluctuation in 
a dark energy component. In the context of a "vanilla" LCDM background cosmology, a non-zero 
zu signals a departure from GR or a fluctuating cosmological constant. Using a phenomenological 
model for the time evolution zu = zuoPoe/ Pm which depends on the ratio of energy density in the 
cosmological constant to the matter density at each epoch, it is shown that the observed cosmic 
microwave background (CMB) temperature anisotropics limit the overall normalization constant to 
be —0.4 < zuo < 0.1 at the 95% confidence level. Existing measurements of the cross-correlations of 
the CMB with large-scale structure further limit zuo > —0.2 at the 95% CL. In the future, integrated 
Sachs-Wolfe and weak lensing measurements can more tightly constrain zuo, providing a valuable 
clue to the nature of dark energy and the validity of GR. 



I. INTRODUCTION 



There is a gap in our understanding of the cos- 
mos: we do not understand the cause of the acceler- 
ated cosmic expansion and the source of most of the 
energy of the universe. In the context of Einstein's 
General Relativity (GR), the evidence suggests the ex- 
istence of some form of dark energy, such as a cos- 
mological constant or quintessence. (For a review see 
Refs. p], 0, H[ and references therein.) Alternatively, 
novel gravitational phenomena on the largest, cosmolog- 
ical scales may be res pon sible for the observations (e.g. 
Refs. 0, H, i, 0, !, S M El)- One of the outstanding 
challenges to physics and astronomy is to discover the 
answer to this problem. Our goal is to reformulate the 
degeneracy between dark energy and novel gravitational 
phenomena into a simple formalism. In doing so, we aim 
to extend the parametriz ed p ost-Newtonian formalism 
(see also Refs. [ll [H [3 El Ell), which has served so 
well to test GR within the solar system, to cosmological 
scales. 

The plan of the paper is as follows. In section II we mo- 
tivate and introduce the new cosmological parameter, zu, 
making a connection with the solar system PPN parame- 
ters. We show that this parameter may characterize new 
gravitational physics, but is also degenerate with shear 
perturbations which we attribute to the dark energy. In 
section III we propose a model for the parameter, zu, de- 
tail our procedure for the computation of cosmological 



observables, and present a sample of the consequences of 
zu ^ 0. Our procedure is not unique but has practical ad- 
vantages. In sections IV and V we present constraints on 
zu and forecasts using the cosmic microwave background 
(CMB), cross-correlations with large-scale structure, and 
weak lensing. A summary with conclusions is given in 
section VI. We include an appendix explaining how dark 
energy with w = — 1 can support fluctuations. 



II. COSMIC PPN 

General Relativity (GR) tells us precisely how mass 
and energy curve spacetime. In other theories of gravity, 
the relationship between between matter and curvature 
is different. By presuming GR and looking for contradic- 
tions we can hope to test the consistency of GR. Post- 
Newtonian variables are used to quantify the behavior 
of gravity and departures from General Relativity. (See 
Refs. pj, [3 for a thorough review.) The Eddington- 
Robertson-Schiff metric 
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can be used to parametrize the way a point mass, M, 
curves space. The Newtonian potential sets Newton's 
constant so that it is convention to set a = 1 at the 
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present day. A more detailed metric and system of 
parameters is possible, allowing for the motion of the 
source, though is unnecessary at this stage. The param- 
eter (3 measures the amount of non-linearity in gravita- 
tional superposition, and 7 measures how much curvature 
is produced per unit rest mass. In GR, f3 = 7 = 1. The 
most recent limits are (3—1 = 1.2 (±1.1) x 10 -4 from 
lunar laser ranging [19], and 7 — 1 = 2.1 (±2.3) x 10~ D 
measured from the time-delay of signals from the Cassini 
spacecraft [2(|. Einstein's General Relativity is in excel- 
lent agreement with observations within the solar system. 

These measurements place very tight constraints on 
extensions of GR. However, it is possible to imagine that 
the underlying theory of gravitation resembles GR on 
solar system scales, but departs on larger, cosmological 
scales. The post-Newtonian variable 7 is very close to 1 
within the solar system, but may depart on larger scales 
or vary with cosmic age. With the notion in mind that 
dark energy phenomena may be a result of a departure 
from GR, it is tempting to speculate that a departure 
from GR hinges on the local ratio of dark energy to the 
matter density. Using standard cosmological numbers, 
this ratio is 

(Pde)/( P m) = n DE ^l/(0.1M G /pc 3 ) « 10- 6 . (2) 

07TG 

Consequently, one may speculate that an improvement in 
solar system tests by an order of magnitude may reveal 
faults in GR. Alternatively, it may be easier to search 
for such faults using cosmological phenomena, such as 
the CMB, the growth of large-scale structure d,[ll|, and 
gravitational lensing [U |22t |23| . 

To adapt the post-Newtonian formalism to cosmology, 
we replace the Eddington-Robertson-Schiff metric with 
the perturbed FRW metric 

ds 2 = a 2 (r) [-(1 + 2^)dr 2 + (1 - 2<j>)dx 2 ] . (3) 

Hereafter we employ the notation and conventions of Ma 
& Bertschinger (MB) (24 |. It is convenient to work in 
the conformal-Newtonian gauge if we wish to maintain 
the simple connection to the Newtonian potential as ap- 
pears in the Poisson equation. If the gravitating source 
is a point mass, then we mean to identify ip — —aGM jr 
and <j) — —jGM/r. so that measurements of the New- 
tonian and longitudinal gravitational potentials can be 
interpreted as a test of the strength of gravity (a) and 
the amount of curvature produced per unit rest mass (7). 
Hence, we introduce a new parameter, zu 1 implicitly de- 
fined by the equation 

^ = (1+^, (4) 

which describes the departure from GR in the cosmolog- 
ical context. In terms of PPN parameters, vj = ^ — 1 
which translates to w ~ 1 — 7 for a weak departure from 
GR. Reexpressing the results of Ref. [23| on the com- 
parison between lensing and dynamical measurements of 
galaxy masses, then w — 0.02 ± 0.07 (68% CL) on tens 



of kiloparsec scales. The same analysis can be extended 
to a few hundred kiloparsec scales by making use of both 
strong and weak lensing measurements towards galaxy 
clusters combined with dynamical mass measurements 
of same clusters. Using the comparison between weak 
lensing and X-ray based masses in Ref. [25|, we find 
zu = 0.03 ± 0.10 (68% CL). At megaparsec scales cor- 
responding to cosmological observations, unfortunately, 
there are no useful limits placed on zu. 

A difference in the two gravitational potentials can also 
occur in GR when the gravitating source has anisotropic 
stress or shear. Specifically, the off-diagonal components 
of the perturbed Einstein equation yield 

k 2 {ct) -ip) = l2nGa 2 (p + p)a, (5) 

where p, p are the energy density and pressure of the fluid 
component giving rise to the shear. Of course, there are 
other sources of shear, due to the contribution of rela- 
tivistic species such as photons or neutrinos. Including 
the effects due to a departure from GR we obtain 

k 2 {cj)-tl)) = \2TiGa 2 {p+p)a\ 1 ^ -zuk 2 cj). (6) 

Note that the difference between the two potentials, as 
well as individual shear components, are gauge-invariant. 
In the synchronous gauge, the potentials are 

(p — rj — Ha, t/j — a + Ha, a — —r»(h + 6rj) (7) 

2£r 

where the dot indicates a derivative with respect to con- 
formal time. Consequently, the shear parametrization 
equation may be recast in terms of the evolution equa- 
tion for a: 

a + (2 + zu)Ha - (I + zu)r] = -12irGa 2 (p + p)<r| 7>1/ . (8) 

Hence, the static PPN-like relation between potentials 
in the conformal-Newtonial gauge becomes dynamical in 
the synchronous gauge. In the concordance cosmological 
model, which includes Einstein's cosmological constant 
within GR, photon and neutrino shear are negligible in 
the late universe. Hence, an observed difference between 
the potentials signals some non-standard behavior. 

This behavior may be seen in gravitational phenomena 
such as the deflection of light. The deflection angle expe- 
rienced by a beam of light moving through a potential <fi 
now takes the form 89 = —2{l+w)4>. In the cosmological 
setting, we can find the lensing potential for the deflec- 
tion of photons by adapting the results of Acquaviva et 
al 26], which amounts to replacing their S with zacf), to 
obtain 

<p(K)= d X g'(x)[C2 + ^)<f>(n,x)}- (9) 
Jo 

Here x is the comoving distance from the observer, 
g'ix) = xf*°° dx'(l - x/x')n( X '), and n(x) gives the 
distribution of background sources. This means the pat- 
tern of lensing of background galaxies, CMB, or any field 
of photons, will be influenced by vu ^ 0. 
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Similarly, the integrated Sachs- Wolfe effect (ISW) is 
affected by w. 

- ¥ {n)\ I sw = J d X ((2 + ™)</>) tT (n, X )- (10) 

Hence, combined measurements of the CMB with probes 
of large-scale structure, which are sensitive to cj> through 
the Poisson equation, can be used to isolate the effects of 

VJ. 

Our focus on the gravitational slip between the po- 
tentials is motivated in part by the recent work of 
Bertschinger |27j, as well as numerous investigations of 
scalar-tensor and gen eraliz ed g ravitational theories {e.g. 
[H, [H, [H ml j}! if) and dark energy with 
anisotropic stress [3J, l35l ] . Our intention is to use these 
existing tools and techniques, not to explore exotic theo- 
ries and models, but to examine simple departures from 
GR with a cosmological constant. 



III. COSMIC MODEL 

It is evident that w can be used to parametrize de- 
partures from Einstein's GR. A specific prediction for 
the gravitational slip in the cosmological context is ob- 
tained for a given theory. For example, the difference 
between the two conformal-Ncwtonian gauge potentials 
is 4> — ip = S(liiL,fj) in a theory of gravity described 
by a generalized Lagrangian L(R, ip) [2a |. where R is the 
Ricci scalar curvature and ip is the gravity-coupled scalar 
field. The relation between a fourth-order gravitational 
theory and PPN parameters 7, /3 is given in Ref. [36| . 
The gravitational slip in a scalar-tensor theory includ- 
ing a scalar coupling to the Gauss-Bonnet Lagrangian 
is analyzed in Ref. 37]. In the Dvali-Gabadaze-Porrati 
(DGP) brane- world model of g ravity [38|,[39[ , recent stud- 
ies of structure formation [401 ] derive a gravitational slip 
(j)-nl) = 8irGa 2 p(5 m + iHV m / k 2 ) / (ik 2 . Here (5 is a di- 
mensionless coefficient which relates the five-dimensional 
cross-over scale, r c , to the present-day Hubble constant, 
whereby w « 0(1/ (3). More examples can certainly be 
found. We would like to find a phenomenological descrip- 
tion which captures the basic effects. The advantage we 
see is that we can streamline the search for non-GR be- 
havior; if we find evidence for w ^ then we can focus 
on particular theories. 

In general, we expect dark energy effects and the cor- 
responding gravitational slip to turn on at late times. 
Hence, we propose a simple extension of the concordance, 
ACDM cosmological model. 1) We assume that the 
dark energy phenomena resembles a cosmological con- 
stant such that the background evolution of the expan- 
sion scale factor is driven by a component with equation- 
of-state w = —1. This maintains consistency with all 
observational constraints based on classical distances. 2) 
We assume an evolution for the gravitational slip which 



grows large as the role of dark energy grows: 
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Our immediate goal is to constrain vjq. Based on our es- 
timate of the local ratio between dark energy and matter 
we expect \wo\ < 10 in order to be consistent with 
the solar system limits on 7 [2(j. Next, using the approx- 
imation that Pde/pm ~ 10 -4 on the scale of a galaxy 
halo, then the limits of Refs. [H SH] give \w \ < 100. 
The cosmological ratio, however, only grows significant 
at late times as the inferred dark energy comes to domi- 
nate. The source of zu, the slip between <fi and ip, may be 
interpreted as due to a departure from Einstein gravity 
or, alternatively, due to a dark energy component with 
anisotropic shear. 3) We employ a subset of Einstein's 
equations to evolve linear perturbations under the pre- 
sumed non-GR theory, to be defined shortly. Hence, we 
are making the implicit assumption that a portion of the 
new field equations are equivalent to GR with an addi- 
tional shear. 

We pause to emphasize the distinction between dark 
energy with shear and departures from GR. Although 
the slip between the gravitational potentials can arise 
from either source, non-GR theories of gravity are not 
generally equivalent to GR and dark energy with shear. 
Our implementation of w represents a only portion of 
the non-standard behavior resulting from different gravi- 
tational theories. Since observations already tightly con- 
strain gross departures from GR, we do not expect to 
be missing much with our parametrization. To keep the 
phenomenology simple, we restrict attention to w = — 1. 
(As well, dark energy or modified gravity with w ^ — 1 
would surely be detected through classical tests of cos- 
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FIG. 1: The CMB anisotropy spectra under recipe Rl are 
shown for different values of zuo- R3 produces identical spec- 
tra. The backgrounds for all models are identical: fi m = 0.35, 
h = 0.65 and dark energy with w = — 1. 
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mology before the effects of shear are manifest.) 

Different theories of gravitation predict different re- 
lations between stress-energy and spacetime curvature. 
Although we presume a A-like background evolution, the 
equations guiding cosmological perturbations must differ. 
In the absence of a specific theory of gravity, we intend to 
use the following recipe (Rl): 3.i) Evolve radiation and 
matter fluid sources using conservation of the individual 
stress-energy tensors, as in the standard case; 3.ii) Evolve 
w according to (TTj) ; 3.iii) Enforce the slip between <j> and 
ip by evolving a according to ([8]); 3.iv) Evolve the time- 
space perturbed Einstein equation k 2 f\ = A-nGa 2 {p + p)9 
where 9 is due to the radiation and matter sources only, 
not dark energy. 

Recipe Rl describes a theory of gravity in which the 
time-space equation is unchanged, but the gravitational 
slip between <j> and ip evolves as given by w. Alterna- 
tively, this describes a dark energy with shear perturba- 
tions but no fluctuations of the momentum density. Im- 
plicitly, there must be energy density and pressure per- 
turbations of the dark energy. (Note that CMBfast 
enforces energy conservation by applying the time-time 
perturbed Einstein equation, h = 2(k 2 ij + 4TrGa 2 Sp)/H, 
at each step. Instead, we use the definition of a to get 
h.) 

We have also considered alternate recipes. For recipe 
R2, the time-time perturbed Einstein equation is evolved 
in place of the time-space equation. This corresponds to 
a theory of gravity in which the time-time equation is 
unchanged; alternatively, the dark energy has no energy 
density fluctuations, but does have momentum density 
and pressure perturbations. The derivative fj is obtained 
from the definition of a. Finally, for recipe R3, we drop 
the evolution of a and let the synchronous gauge vari- 
ables h and rj evolve as in the standard case. However, 
for matters of calculating the CMB anisotropics we con- 
vert to the conformal-Newtonian gauge with <p = rj — Tia 
but then artificially set ip = (1 + m)(f>. 

We have modified a version of CMBfast [42j to evolve 
perturbations according to these schemes. The resulting 
CMB anisotropy spectra are shown in the following fig- 
ures. All the models are spatially flat with f2 m = 0.35, 
h = 0.65, and normalized to have the same small angular- 
scale power spectra as the corresponding ACDM model. 
In Figure [1] we show the anisotropy power spectrum for 
different values of vjq under recipe Rl. The alternative 
model, R3, produces identical spectra. Empirically, we 
see that for small, negative values of wq close to —0.1 
there is a cancellation between the ISW and SW ef- 
fects which leads to a slight reduction in the amplitude 
of the lowest multipole moments. For these cosmologi- 
cal parameters, the minimum value of the quadrupole is 
achieved with too ~ —0.15. For large values of | vjq +0.15 1, 
the reduction turns into an enhancement. In Figure [5] we 
show the anisotropy power spectrum with wq = —0.1 
for different recipes. All three recipes produce similar 
results. 

The growth of structure is also modified by the gravita- 
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FIG. 2: The CMB anisotropy spectra under different recipes 
are shown for voq = —0.1. 



tional slip. According to recipe Rl, the cold dark matter 
density contrast S c evolves as 

8 C + H6 C = k 2 {l + m){Ha-rj). (12) 

where we have used 5 C = —h/2, equation ([8|), and we 
have assumed that the shear due to radiation and neutri- 
nos is negligible, as occurs beginning late in the matter- 
dominated era. Recipe Rl maintains the validity of 
the time-space perturbed Einstein equation, whereas the 
time-time equation is implicitly modified. If the mod- 
ifications of the time-time equation are absorbed into a 
dark energy density perturbation, whereby k 2 rj — Hh/2 = 
—AnGa 2 (5p+8pDE), then the CDM density contrast evo- 
lution equation becomes 

S c + H5 C = ATtGa 2 (5 P + 5p DE ){l + vj). (13) 

The effective dark energy density perturbation on the 
right-hand side is familiar from studies of quintessence 
dark energy. The novel feature of this equation is the 
1 + vj factor which can lead to faster (slower) growth for 
vo > (< 0). (Also see Refs. 0, 0, g^.) 

The modification of standard cosmological perturba- 
tion theory to incorporate the effects of gravitational slip 
is not unique, as described above. There may be recipes 
other than Rl-3 to include w due to a departure from 
GR. Any particular non-GR theory makes a specific pre- 
diction for the behavior of w as well as other effects. In 
the appendix, we study the behavior of a dark energy 
component with w ~ — 1 which also generates a gravita- 
tional slip. However, as we have shown, the differences 
in the three recipes are small when wq is small. Thus, 
despite the degeneracy, we find a practical advantage in 
this simple modeling of departures from GR. 
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IV. CMB ANISOTROPY SPECTRUM 

Let us now investigate to what extent zuq can be con- 
strained from current CMB observations. As shown in 
the previous section, a change in zuq will mostly affect the 
large angular scale CMB anisotropies through the ISW 
effect, leaving the acoustic peak structure unchanged. 
For a first study, we restrict the analysis to a flat, "con- 
cordance," ACDM model with Q A = 0.76, Q b = 0.02, 
h = 0.73. We also make the choice of adiabatic and 
scalar inflationary perturbations with a spectral index of 
n s = 0.958 and we fix the optical depth to r = 0.08 
as suggested by Wilkinson Microwave Anisotropy Probe 
(WMAP) 3-year anisotropy and polarization measure- 
ments. We then let wq vary and compare the CMB 
power spectra with the WMAP three year temperature 
anisotropy data 0, S3, M, HI (WMAP3). The likeli- 
hood is determined using the October 2006 version of the 
WMAP likelihood code available at the LAMBDA website 
f50j ] . In Figure [3] we plot the likelihood for tu (under 
recipe Rl) from WMAP. As we can see, the likelihood 
peaks at vdq ~ —0.1 although -cuq = is consistent with 
the data: the 95% CL range is -0.28 < ro < 0.05. The 
preference for negative values of vjq is a direct conse- 
quence of the low measured CMB quadrupole. Hovewer, 
as shown in the previous section, very large negative val- 
ues such that \w§\ > 0.5 enhance the quadrupole, making 
the model dramatically incompatible with the data. 

It is also interesting to investigate possible variations 
of zuq by relaxing the assumption of w = —1 dark en- 
ergy component. Let us indeed substitute the cosmolog- 
ical constant with a perturbed fluid with constant equa- 
tion of state p = wp and sound speed c 2 s . An equation 
of state w 7^ — I and perturbations in the fluid would 
also change the large angular scale anisotropy. As it is 



well known (see e.g. [5l|,[52j]) a geometrical degeneracy 
makes virtually impossible any determination of w from 
the position of the acoustic peaks in the CMB anisotropy 
spectrum. However, combined measurements constrain 
the equation of state parameter to a conservative range 
-1.2 < w < -0.8 (see e.g. 46]). In Figure [3] we show 
the likelihood for vjq when a dark energy fluid with equa- 
tion of state w = —0.8 is considered. The other param- 
eters have been chosen to reproduce the acoustic peak 
structure of the concordance model through the geomet- 
rical degeneracy relation. We therefore fix Q m = 0.298, 
h = 0.655 and flf, = 0.0495. We consider two possi- 
bilities for the sound speed parameter = 1, 0. As 
we can see, the likelihood for wq is only moderately af- 
fected: we obtain —0.38 < wq < 0.05 for <? a = 1 and 
-0.35 < vj q < 0.09 for c 2 s = at 95% CL. We can 
therefore claim that current observations provide a con- 
servative bound of -0.4 < w < 0.1 at 95% CL. This 
constraint is stronger than the PPN, lensing, and x-ray 
constraints on zuo, based on our model (|II[) . because 
Pdm "C Pm on megaparsec scales and below. 
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FIG. 4: Likelihood distributions for zu due to WMAP 3- 
year CMB temperature measurements within the framework 
of different dark energy models. The curve with dots is 
for ACDM, as in the previous figure. The solid and dotted 
curves are for a dark energy model with w = —0.8, with pa- 
rameters chosen to reproduce the acoustic peak structure of 
the ACDM model, but with dark energy perturbation sound 
speed <? a — 1,0. The sound speed does not have a strong 
effect on the likelihood range. Allowing for a range in equa- 
tion of state —1.2 < w < —0.8, then we obtain a conservative 
bound -0.4 < uj < 0.1 at 95% CL. 



FIG. 3: Likelihood distribution for ujq due to WMAP 3-year 
CMB temperature measurements for a flat, "concordance," 
ACDM model with Q A = 0.76, Q b = 0.02, h = 0.73, adiabatic 
and scalar inflationary pertubations with a spectral index of 
n s = 0.958, and an optical depth of r = 0.08. The 95% CL 
range is —0.28 < vjq < 0.05 with a peak at a slightly negative 
value of vjq owing to the decrease in large-scale power. 



EXPECTATIONS FOR ISW AND WEAK 
LENSING MEASUREMENTS 



Post-GR gravitational slip would also leave an imprint 
on the cross-correlation between the CMB and large- 
scale structure. We define the two-point angular cross- 
correlation between the temperature ISW anisotropy and 
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the dark matter fluctuation as (see e.g. |53|): 

C x {0) = (A ISW (ix)6 LS s(i2)), (14) 

where the angular brackets denote the average over the 
ensemble and 9 = |t\ — 72 1. For computational purposes 
it is convenient to decompose C x (9) into a Legendre se- 
ries such that, 

C x (f)=£^Q*fl(coB(0), (15) 
1=2 

where Pi(cos9) are the Legendre polynomials and C x is 
the cross-correlation power spectrum given by [54l . [Eg : 

C*=4ir^JfAlir W (k)lt SS (k), (16) 

where is the primordial power spectrum. The inte- 
grand functions tf sw (k) and lf" ss (k) are defined respec- 
tively as: 

It ss (k) = b J <S>{z)8 k m {z)3i[kr{z)]dz, (18) 

where 4>k and 5^ are the Fourier components of the grav- 
itational potential and matter perturbation, respectively; 
<£> is the galaxy survey selection function; ji[kr{z)] are the 
spherical Bessel functions; r(z) is the comoving distance 
at redshift z and k(t) = k(t)g?t is the total optical 
depth at time r. 
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FIG. 5: Contributions to the integrated Sachs- Wolfe effect as 
a function of the redshift z for different values of xuq. 

A change in vj could therefore change not only the 
value of lj sw but also its sign respect to l[" ss . In Fig- 
ure we plot the values of the integrand e~ K d((2 + 
w)4>k)/dT at a scale of k ~ 0.0016 Mpc -1 as a function of 
different values of zu . As we can see, while the early ISW 
at z ~ 1000 is left unchanged, a value of U7 ~ —0.2 pro- 
duces a change in the sign of the late ISW contribution. 
This could indeed produce a negative cross-correlation 
with galaxy surveys. 



Direct measurements of the cross-power spectrum C x 
are more robust for likelihood parameter estimation since 
these data would be less correlated than measurements of 
C x (9) . We therefore compute the cross-power spectrum 
C x for different values of assuming a galaxy survey 
with a selection function as 

$(z) - z 2 exp[-(z/z) 1 ' 5 ] (19) 

where z, the median redshift of the survey, is z = 0.25. In 
the past years, the WMAP temperature anisotropy maps 
have been cross-correlated with several surveys of Large 
Scale Structure (LSS) distributions and a positive corre- 
lation signal has been detected [H, [UIH, [5t| [g3,[6l|. As 
we can see from Figure |H1 values of wq < — 0.2 would re- 
sult in anti-correlation and a negative angular spectrum 
which is in disagreement with current observations at the 
~ 2cr level. 
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FIG. 6: Cross-correlation angular power spectrum between 
CMB temperature (ISW) and galaxy distribution as a func- 
tion of vjq. A value of ujq < —0.2 changes the sign of the 
ISW, resulting in a negative cross-correlation, in conflict with 
observations. 

In addition to the ISW effect, weak lensing measure- 
ments of the large-scale structure capture modifications 
imposed by vj. Using equation ©, we can write the 
angular power spectrum of lensing convergence k as 

C ? = IJ k2dk t /K ( fc )] 2 P <^ fc ) ' ( 2 °) 

where 

I K (k) = J d X g( X )(2 + ™)ji{kx) ■ (21) 

In Fig. [3 we show the lensing convergence for our stan- 
dard LCDM model and a model in which wq = — 0.1 and 
wq = 0.1 in Eq. (11) assuming that all sources are at a 
fixed redshift of z s = 1. As shown, there is an overall 
change in the amplitude of lensing convergence fluctu- 
ations such that when -cjq = —0.1 (0.1) the fluctuation 
power specrum is lower (higher) by about 8%. To see the 
modification associated with w, the overall normalization 
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FIG. 7: Weak lensing convergence power spectrum for LCDM 
and for several values of O7o assuming that all sources are at 

Z S = 1. 



of the convergence power spectrum, such as through eg, 
must be known to an accuracy better than 8%. 

Detection of modifications from w with lensing statis- 
tics alone is impossible in the same manner that lensing 
measurements of individual foreground objects cannot be 
used to separate cj) from i/j as the measurements are sen- 
sitive to the combination only. To constrain our post- 
GR parameter, as in Ref. (23| . one must combine lensing 
measurements with an independent measurement which 
is sensitive to a different combination of the two parame- 
ters or to just one of the parameters. In the case of galaxy 
lensing, the comparison is made to dynamical mass es- 
timates which are sensitive to 4>. In this context, in the 
case of large-scale structure weak lensing, the weak lens- 
ing statistics must be compared with an additional probe 
of large-scale structure, such as the potentials probed by 
the galaxy distribution. While we have not carried out 
such an analysis here due to the as-yet large uncertain- 
ties in weak lensing measurements, in the future it may 
be possible to further constrain through a combina- 
tion of lensing measurements and galaxy clustering power 
spectrum. 



VI. CONCLUSIONS 

We have introduced a new cosmological variable to 
characterize the degree of departure from Einstein's Gen- 
eral Relativity. The new parameter, w, is the cosmolog- 
ical analog of the parametrized post-Newtonian variable 
7, which measures the amount of spacetime curvature 
per unit mass. Our study was motivated by various dis- 
cussions in the literature on uses of large-scale structure 
measurements to establish departures from GR on scales 
beyond the solar system (e.g., [27]). In the cosmolog- 
ical context, the parameter w measures the difference 



between the Newtonian and longitudinal potentials in 
response to the same matter sources, as occurs in cer- 
tain scalar-tensor and other modified theories of gravity. 
Equivalently, w measures the scalar shear fluctuation in 
a dark energy component. In the context of a "vanilla" 
LCDM background cosmology, a non-zero w signals a de- 
parture from GR or a fluctuating cosmological constant. 

We have parameterized the time-evolution of zu 
through a simple phenomenological model in which zu = 
zuo(l + z)~ 3 Qy)e/Qm- The limit due to the observed 
cosmic microwave background t emp erature anisotropics 
from the WMAP3 data [M SI H HI is -0.4 < w < 
0.1 at the 95% confidence level. The detection of a pos- 
itive ISW effect [H [H [H, [H, M, Ell further limits 
zuq > —0.2. Within the context of our model for zu(z), 
these constraints are tighter than the solar system time- 
delay and galaxy-scale gravitational lensing tests of GR 
by an order of magnitude. 

We have considered possibilities to improve these con- 
straints and have discussed the potential use of both 
ISW and weak lensing measurements of the large-scale 
structure. In both cases, these measurements must be 
combined with additional probes of large-scale structure, 
such as the galaxy power spectrum, to break degenera- 
cies between w and parameters. While existing data are 
limited given the large uncertainities involved, weak lens- 
ing measurements combined with galaxy clustering may 
be an ideal probe of w. Further evidence that \zu\ <C 1, 
or a detection of w ^ 0, will provide useful clues as to 
the validity of GR and the nature of dark energy. 
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APPENDIX A: W=-l DARK ENERGY 
FLUCTUATIONS 

It is an apparent contradiction that a cosmological con- 
stant can fluctuate. Einstein's cosmological term is re- 
garded as a universal constant, leading to a stress-energy 
= — Ag^/8TrG. The fluctuations must vanish by defi- 
nition, 8T{? = 0. To avoid confusion, we will keep the cos- 
mological constant pristine - anything otherwise would 
be a misnomer. However, if the dark energy is a con- 
stant only after some suitable averaging or expectation 
value is obtained, then it is not clear that its fluctuations 
must also vanish. An example which has been explored 
recently is the ghost condensate (62|, consisting of an 
effective field theory which gives rise to a cosmological 
fluid with equation-of-state w = —1, but also supports 
scalar fluctuations and a novel dispersion relation. The 
full implications of such a theory are beyond the scope 
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of our investigation, but serves as a motivating proof-of- 
principle. 

In the absence of a specific model, we describe the 
dark energy as an imperfect fluid with equation-of-state 
w = p/p for the homogeneous, isotropic background, but 
which also supports anisotropic stress or shear, a. Ex- 
amples of dark energy models with anisotropic stress are 
given in Refs. [34l 1.351 l63l l64j . The equations of motion 
for the linearized perturbations are given by 



= -H(l-3w)9- 



6p/Sp -k 2 5-k*a + k^ 



1 + w 1 + w 
8 = -(1 + w) (0 - 3^) - mf^-ivjs (Al) 

which are merely equations (30a, b) of [24j | . Note that 
we neglect vector and tensor perturbations, which may 
play a role in dark energy models with shear. The 
pressure and energy density perturbations are linked by 
c 2 = Sp/Sp. If we consider that V = (1 + w)0 is the 
relevant variable to describe the momentum density per- 
turbation, and II = (1 + w)a describes the shear, then 
the equations become 

V = -H(l-3w)V + c 2 s k 2 6-k 2 Il+(l+w)k 2 ip 
6 = -V + 3(1 + w)4>- 3H(c 2 -w)S. (A2) 

Next, we set w = — 1 and make the reasonable assump- 
tion that the gravitational potentials <f> and ip remain 
finite but that V, 8p, 5p, II do not automatically vanish. 
The equations become 



V = -AHV + c 2 k 2 S - fc 2 n, 
6 = -V-3H((£ + i)6 



(A3) 



There is no longer any source term representing the in- 
homogeneities of the surrounding medium; the fluctua- 
tions only feel the background, cosmic expansion. Of 
course, these fluctuations still feed in to the metric per- 
turbation equations. Note that the set of equations (|A3[) 
become trivial in the case of a scalar field, for which 
II = (5 = V = 0in the limit w — ► — 1. 

The system of equations (|A3[) can be expressed as a 
single, second order equation for 5: 

8 + (7 + 3c 2 )nS+ [(I2H 2 + m'){\ + c 2 ) + c 2 k 2 ] 5 = k 2 Il. 

(A4) 

In the absence of shear, the short-wavelength modes are 
stable provided c 2 > 0. The long- wavelength modes de- 
cay as a~ 3 , a~ 3 ( 1+c =- ) in a radiation-dominated epoch, 
and as a~ 7 / 2 , a~ 3 ( 1+c =) under matter. Any such fluc- 
tuations are of little interest since the decay is so rapid. 

A wider variety of behavior is possible with the addi- 
tion of shear. The sound speed c s no longer serves as the 
relevant quantity guiding the propagation and stability of 



perturbations. Now, c 2 8 > II is required for stability, and 
the high-frequency phase velocity is v — (c 2 — H/5) 1 / 2 . 
This corresponds to a time and scale-dependent cosmic 
post-GR parameter 

w = -12TrGa 2 pU/k 2 cj). (A5) 
The propscrtics of the shear, and therefore the cosmic 
post-GR parameters, are dictated by the details of any 
particular model. The quartic dispersion relation for the 
excitations in the ghost condensate can be modeled by 



n : 



M 2 



S. 



(A6) 



where M is the energy scale below which the effective 
field theory is valid. Stability of modes within the horizon 
requires c 2 > 0, in which case the perturbations decay 
rapidly. Next, a phenomenological evolution equation 
has been proposed by Hu [64]], whereby 



n + 3HIT= -c 2 



(A7) 



(when translated into our notation) in the conformal- 
Newtonian gauge. The viscosity sound speed c V i S mod- 
ulates the source of the shear, which otherwise decays. 
However, a phenomenological model which reflects equa- 
tion J2]), whereby the shear-effects turn on as the lo- 
cal dark-energy density overtakes the matter, can be 
achieved if 



II = — ZUn5. 

3 



(A8) 



In this case, w ~ vjqSpde/^ptot on scales within the 
horizon, which resembles our model (TIT]) . We can also 
make a connection with Hu's formula by plugging 
into (IA3D. such that 



n + 3W(i 



2 IcHnnJo 
c = - -— II 



1 



3 din i 



-ZUqV. 



(A9) 



This corresponds to a viscosity sound speed c 2 is = 
—wq(1 + w)/8, but a damping rate which is now variable. 
With c 2 < — 1 and a constant or growing parameter vjq, 
the anisotropic scalar shear can actually grow. It is not 
clear what would be an appropriate choice of initial con- 
ditions for fluctuations in this dark energy component, 
since adiabaticity requires 6 = V = 0. Initial conditions 
which are set by some measure of the relative contribu- 
tion to curvature perturbations, 5R = —8ttG(Sp — 38p), 
must be finely tuned in order to force a non-negligible 
gravitational slip at late times. In this investigation, 
however, we have focused directly on the slip rather than 
fluctuating dark energy models. 
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